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ABSTRACT 

We performed a systematic analysis of all INTEGRAL observations from 2003 to 2009 
of 14 Supergiant Fast X— ray Transients (SFXTs), implying a net exposure time of about 
30 Ms. For each source we obtained lightcurves and spectra (3 — 100 keV), discovering sev- 
eral new outbursts. We discuss the X— ray behaviour of SFXTs emergin g from o ur analysis 
in the framework of the clumpy wind accretion mechanism we proposed dPucci e t al. 200S^. 
We discuss the effect of X— ray photoionization on accretion in close binary systems like 
IGR J16479-4514 and IGR J17544-2619. We show that, because of X-ray photoioniza- 
tion, there is a high probability of formation of an accretion disk from capture of angular 
momentum in IGR J16479— 4514, and we suggest that the formation of transient accretion 
disks could be responsible of part of the flaring activity in SFXTs with narrow orbits. We also 
propose an alternative way to explain the origi n of fla res with peculiar shapes observed in 
our analysis applying the model of ,Lamb et al.1 (|1977|) . which is based on the accretion via 
Rayleigh— Taylor instability, and was originally proposed to explain type II bursts. 

Keywords: X-rays: binaries -X-rays: individuals: IGR J16479-4514,XTE J1739-302/IGR 
J17391-3021,AXJ1841.0-0536/IGRJ18410-0535,IGRJ18483-0311,SAXJ1818.6-1703, 
IGR J16418-4532, AX J1820.5-1434, AX J1845.0-0433, IGR J16195^945, IGR J16207- 
5 1 29, IGR J 1 6465-4507, IGR J 1 7407-2808, XTE J 1 743-363 . 



1 INTRODUCTION 

Supergiant Fast X-ray Transients (SFXTs) are a sub-class of High 
Mass X-ray Binaries (HMXBs) discovered by the INTEGRAL 
satellit e in the last seven years, during the Galactic plane moni- 
toring jSguera et alj 2005). SFXTs host an OB supergiant and an 
accreting compact object, and show a sporadically X— ray tran- 
sient emission composed by many flares reaching a luminosity of 
10»e _ 10=*^ erg s~ \ with flare durations of ~ 10^ - lO" s. For most 
of their lifetime, SFXTs accrete at an intermediate level, showing 
an X— ray luminosity of 10^^ — lO'^^ erg s~^, as di s covered by the 
Swift fXSH monitoring (see e.g. ISidohetalJl2008l ; iRomano et al.l 
l2009t) . The quiescent level (a bout 10^^ erg s~ ^ ) has been observed 
only in a few so urces (see lin't Zand l2005l : iLevder et al.l 1 20071 : 
IBozzo etal.ll2010t) . Thus SFXTs show a large dynamic range of 
about 10^ — 10^. Many accretion mechanisms have been proposed 
to explain the transient behaviour, lin't Zandl j200^ suggested that 
the large dynamic range could be p roduced by t he acc retion of 
dense clumps from the donor wind. ISidoli et^ ( l2007h invoked 
the presence of an equatorial wind component, denser than the po- 



lar wind component and inclined with respect to the orbital plane 
of the compact object to explain the transient periodi c emission 
of IGR J11215-5952. IBozzo. Falanga & Stellj ( 12008 j) proposed 
that SFXTs host a magnetar with large spin period (~ 10'^ s): the 
changes in X— ray luminosity are ascribed to the gated accretion 
mechanisms, where changes in the reciprocal positions of accre- 
tion, magnetospheric and corotation radii lead to different accretion 
regimes. 

In this paper we report the results from the analysis of INTE- 
GRAL data of 14 SFXTs, for a total exposure time of ~ 30 Ms. 
The results obtained here are dis cussed considerin g the structure 
of the clumpy supergiant winds l lDucci et"aill2009l) . the effect of 
X— ray photoionization of the outflowing wind b y the compact ob- 
ject (in the framework of the lHo & Arons|[T987l accretion model), 
and the pos sible formation of transient ac cretion disks, as those 
proposed bv lTaam. Brown & Frvxelll ( Il988l) to reproduce the flares 
from EXO 2030+375. 



2 OBSERVATIONS AND DATA ANALYSIS 



* E-mail: lorenzo@iasf-milano.inaf.it 



The INTEGRAL observatory, launched in October 2002, carries 
3 co-aligned coded mask telescopes: the imager IBIS (Imager on 
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Figure 1. IBIS/ISGRI lightcurves of IGR J16479-4514, XTE J1739-302, IGR J17544-2619, AX J1841.0-0536 (20 - 40 keV). A binned time correspond- 
ing to the ScW duration (~ 2000 s) has been used. 



Table 1. Total exposure time Toxp, duration and number of ScWs where 
the source is detected with a significance greater than 5. 



Source Toxp (d) Duration (d) # ScWs # Outbursts 

S/N > 5 S/N > 5 

38 
18 
14 
4 
13 
11 
23 
2 
7 
3 
9 
2 
3 
7 



Board the INTEGRAL satellite: lub"ertini et alj2003h . sensitive from 
15 keV to 10 MeV, the spectrometer SPI (SP ectrometer on INTE- 
GRAL, 20 keV-8 MeVi lVedrenne et al.ll2003l) . and the two X-ray 



monitors JEM-X l and JEM-X2 (Joint European X— ray Monitor; 
iLund et aljr2003l) . sensitive in the energy range 3—35 keV. IBIS 
is composed of the low-energy detect or ISGRI (INTEGR AL Soft 
Gamma Ray Instrument; 15-600 keV; lLebmn et alj|2003l) and the 
Csl layer PICsIT ( Pixellate d Imaging Caesium Iodide Telescope; 
175 keV-10 MeV; iLabanti et al. 2003) . /ATEGiML observations 
are divided in pointings called Science Windows (ScWs), which 
have a typical exposure of 2 ks. 

We analysed all the public and our private data, between 2003 
and 2009, where the SFXTs IGR J16479-4514, XTE J1739-302, 
AX J1841.0-0536 and IGR J17544-2619 were within 15° from 
the center of the field of view. This resulted in 14426 Science Win- 
dows, corresponding to a total exposure time of ~ 30 Ms. We 
considered in our analysis also the other SFXTs (and candidate 
SFXTs) observed in the ScWs selected (see Table[T). We analysed 
IBIS/ISGRI and JEM-X data using the Off-line Scientific Analysis 
package OSA 8.0 (Goldwurm et al. 2003). For the spectral analy- 
sis, which was performed with XSPEC (ver. 11.3), we added a 2% 
systematic error to both IBIS/ISGRI and JEM-X data sets. 

For each source reported in Table [T] we extracted the 
IBIS/ISGRI lightcurve at a ScW time resolution (~ 2000 s), 
in the energy range 20—40 keV, and we considered only the 
pointings where the sources are detected with at least a 5a sig- 
nificance. For each source. Table [T] reports the number of out- 
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Figure 2. Joint JEM-X plus IBIS/ISGRI counts spectra, together with residuals in units of standard deviations (left panels) and unfolded spectra (right panels) 
for IGR J16479-4514. XTE J1739-302. and IGR J17544-2619 (see Tableg). 



burstfl we identified. For each outburst observed, we also extracted 
IBIS/ISGRI lightcurves (18-60 keV) with a bin time of 50 s, and 
an IBIS/ISGRI spectrum in the energy range 18—100 keV. We 
found that the outbursts are characterized by a flaring activity, typ- 
ical of SFXTs. Due to its smaller field of view, we were able to 
extract the JEM-X spectra in only a few pointings. 



^ We defined outburst as an X— ray emission detected by IBIS/ISGRI with 
significance > 5cr, and separated by adjacent outbursts by at least ~1 day 
of inactivity (source below the IBIS/ISGRI threshold of detectability). An 
outburst can be composed by one or a series of flares with a typical duration 
of ~ 10^ - 10-* s each. 



3 RESULTS 

With IBIS/ISGRI we discovered several previously unno- 
ticed outbursts from XTE J1739-302, IGR 116479-4514, 
IGR J17544-2619, IGR 118483-0311, AX J1845.0-0433, 
IGR 116195-4945, IGR J16465-4507, IGR J17407-2808, 
IGRJ16207-5129, and IGR J16418-4532. We reported these out- 
bursts in Table |2] and [3] together with the mean flux and best fit 
parameters of the IBIS/ISGRI spectra. 

Figure [T] displays the 20-40 keV IBIS/ISGRI lightcurves 
of IGR 116479-4514, XTE J1739-302, IGR 117544-2619 and 
AX J1841.0-0536 collected from 2003 to 2009, where each "de- 
tection" refers to the average flux observed during each ScW, when 
the sources are detected at a significance > 5cr. Solid boxes rep- 
resent the observations where the sources are not detected (< 5a). 
For most of the time each source is not significantly detected. 
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Table 2. Summary of the new flares discovered of 7 confirmed SFXTs (IBIS/ISGRI data). We fitted the spectra with a powerlaw or a bremsstrahlung model. 



Source 


Date 






Mean Flux 


kT 


(kevyr 


xl (d-o.f ) 




(UTC) 






(erg cm ^ 


s-i) 








XTEJ1739-302 


2009 Apr. 


8, 16:49- 


18:46 


4 X 10" 


10 


r = 


9 a +0.5 


0.99 (10) 


IGR J 16479-45 14 


2003 Feb. 


2 19:07-22:30 


2.5 X 10 


-10 


r = 


2 g+0.5 


0.76 (9) 




2004 Aug 


9, 02:26- 


04:24 


3.0 X 10 


-10 


kT -- 


= 41 

-17 


0.52 (13) 




2004 Aug 


20, 07:27 


-12:58 


2.6 X 10 


-10 


r = 


— 0.5 


0.73 (10) 




2004 Sep. 


10, 01:14- 


-02:12 


2.5 X 10 


-10 


kT - 


= 19+" 


0.4(12) 




2009 Mar. 


1, 16:56- 


17:54 


10-3 




r = 


9 7+0.7 
-0.6 


1.2(9) 




2009 Apr. 


6, 6:16-13:09 


7.2 X 10 


-10 


kT 


= 34+^ 


1.45 (12) 



IGR J 17544-26 19 



2004 Feb. 27, 14:18-14:45 
2006 Sep. 20,09:58-13:47 



4.2 X 10-1" 
6.9 X 10-10 



r = 

kT : 



3.1 



18 



+0.9 
-0.8 
+29 



1.27(8) 
1.38 (7) 



IGR J 18483-03 11 2004 Nov. 3, 00:05- 


00:31 


3.3 X 10- 


10 


kT 


_ 99+48 
— _ii 


1.01 (11) 


2005 Oct. 16, 07:34- 


-08:30 


1.7 X 10- 


10 


kT 


_ 91 +39 
- ^1 -10 


1.19 (11) 


2005 Oct. 19, 20:48- 


-20, 18:40 


3.1 X 10- 


10 


r = 


. 9 O+0.3 


1.42 (10) 


2006 Apr. 25, 14:04 


-14:33 


5.0 X 10- 


10 


kT 


= 40+^1 
-21 


1.54(11) 



AX J1845.0-0433 


2006 Sep. 3, 17:52- 


18:28 


4.6 X 10- 


10 


r 


- 2 4+"'^ 

— ^-^-0.5 


0.83 (11) 


IGRJ16195-4945 


2004 Aug. 20, 04:40 

2005 Feb. 18, 13:25- 


-06:27 
-14:22 


2.6 X 10- 

2.6 X 10- 


10 
10 


r 
r 


- 2 0+" '^ 

- 2 0+^'^ 


0.73 (10) 
0.61 (8) 


IGRJ16465-4507 


2004 Aug. 9, 22:12- 


10, 02:53 


2.3 X 10- 


10 


r 


_ 9 c +0.5 


0.86 (10) 



AX J1841. 0-0536 




Energy (keV) 

Figure 3. Average IBIS/ISGRI counts specfi-a of AX J1841.0-0536 (see 
Table|4j, and residuals (lowest panels) in units of standard deviations. 
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No evidence of spectral variability has been found between 
different outbursts of these four SFXTs. To achieve the best statis- 
tics, we extracted an average JEM-X-I-IBIS/ISGRI outburst spec- 
trum for each source (Figures[2]and[3}. We did not extract an aver- 
age JEM-X spectrum for AX J1841. 0—0536 because we detected 
the source only in one pointing, with a low flux. The models which 
best fit the average spectra are reported in Table|4l 

3.1 Clumpy wind in IGR J16479-4514 

The results of our analysis of all available INTEGRAL observations 
of SFXTs can be compare d with quantitative expectations from our 
new clumpy wind model toucci et al.ll2009l) . A meaningful com- 
parison can be performed only in SFXTs where a large number of 
flares has been observed. For this reason we will concentrate only 



Figure 4. Folding on a period of 3.3194 d of the ScWs where 
IGR J16479— 4514 is detected with significance > 5cr. The start time of 
the eclipse is to = 54546.742 MJD, and the duration is 0.6 d. 



on IGR J16479— 4514. In order to apply the model, we have to 
establish the stellar parameters of the system and to confirm the 
eclipse duration. 

In Figure |4] we show the IBIS/ISGRI ScWs where 
IGR J 16479— 4514 has been detected in outburst, folded on 
the orbital perio d of the system as suming a zero time to = 
54546.742 MJD ( IBozzo et al]|2008bh . This histogram is consis- 
tent with the presence o f an eclipse with duration At ~ 0.6 d 
jjain. Paul & Duttj2009l) . 

The Roche Lobe radius Rl of IGR J 16479-45 14, adopting 
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Table 3. Summary of the new flares discovered of 3 candidate SFXTs (IBIS/ISGRI data). We fitted the spectra with a powerlaw or a bremsstrahlung model. 



Source Date Mean Flux kT (keV)/r xl (d.o.f) 

(UTC) (ergcm-^s-i) 



IGR J 17407-2808 


2003 Sep. 21,04:27-04:29 


2.6 X 10- 


-9 


kT -- 


— 44 +35 
~ 16 


1.38(10) 




2006 Sep. 20, 12:18-13:47 


4.9 X 10" 


-10 


V = 


1.9+°'^ 


0.55(14) 


IGR J16207-5129 


2004 Jan. 21, 09:23-09:56 


5.8 X 10- 


- 10 


kT = 


= 38.1^7 


1.37(10) 




2005 Feb. 07, 06:48 -07:45 


2.0 X 10- 


- 10 


kT -- 




0.69(11) 




2005 Feb. 11,03:04-04:01 


2.2 X 10" 


-10 


r = 


2 ^ 


0.78 (10) 




2005 Feb. 13, 10:59-11:56 


2.2 X 10- 


10 


r = 


_un s 


1.06(10) 




2005 Feb. 18, 22:41-19, 01:41 


3.2 X 10- 


10 


r = 


^■^-0.3 


1.06(11) 


IGR J16418-4532 


2003 Feb. 3, 12:37-13:13 


4.6 


X 


10- 


-10 


r = 


1.8+;-? 

— 1.1 


1.07 (5) 




2003 Mar. 4, 21:39-22:07 


2.1 


X 


10- 


-10 


kT -- 


= 15+«' 

— 6 


1.12 (14) 




2004 Feb. 18, 03:42-04:16 


3.0 


X 


10- 


-10 


V = 


3.2+^2 

— 1.0 


1.75(14) 




2004 Mar. 20, 21:27-21, 11:24 


2.5 


X 


10- 


-10 


V = 


3.0+°I 

— 0.6 


0.96 (14) 




2004 Aug. 9, 13:57-14:56 


1.9 


X 


10- 


10 


r = 


2.7+!;'^ 

-0.7 


0.70 (14) 




2004 Aug. 24, 13:32-14:56 


2.0 


X 


10" 


-10 


kT -- 


= 20+12 

— 7 


1.00 (14) 




2005 Feb. 13,02:17-15:02 


1.6 


X 


10- 


-10 


kT 


= 15^« 


1.34(14) 




2005 Feb. 22,04:54-06:14 


1.5 


X 


10- 


-10 


kT -. 


OO +23 

- ^^-10 


1.13(14) 




2005 Feb. 24,05:20-07:10 


1.3 


X 


10- 


-10 


V = 




0.87(14) 




2005 Feb. 27, 05:41-Mar. 1, 12:09 


1.9 


X 


10- 


10 


kT 


= 21+5 


1.02(14) 




2005 Mar. 8, 04:35-05:35 


3.2 


X 


10- 


-10 


kT 


= 14+* 


1.16(14) 




2005 Mar. 18, 14:09-14:42 


1.3 


X 


10- 


-10 


V = 


4 2+5-8 


0.62(7) 




2005 Aug. 26, 18:08-27, 08:56 


2.5 


X 


10- 


10 


kT -. 


= 28+lg5 


1.16(14) 




2005 Oct. 2,00:32-01:15 


3.6 


X 


10- 


-10 


V = 


2 6+l'^ 


0.76 (8) 




2006 Aug. 13, 23:48-14, 21:01 


2.9 


X 


10- 


-10 


kT 


= 20+^ 


1.65(14) 




2007 Feb. 19, 03:58-04:39 


3.9 


X 


10- 


-10 


V = 


r, O+0.9 

^■•^-0.8 


0.74 (8) 




2007 Mar. 15, 04:43-05:27 


5.4 


X 


10- 


-10 


kT -. 


= 15+_- 


1.48 (8) 




2007Mar. 30, 02:19-03:01 


3.0 


X 


10- 


-10 


V = 




0.87(5) 




2009 Apr. 7, 10:15-10:43 


4.8 


X 


10- 


-10 


V = 


2.8+n 


0.92(10) 




2009Apr. 7, 21:13-21:41 


4.1 


X 


10- 


-10 


r 9 OS +0.99 
1 - 2.28 _o 3^ 


1.22(9) 



Table 4. Best fit parameters of the average spectra of IGR J16479-4514, XTE J1739-302, and IGR J17544-2619 observed with JEM-X and IBIS/ISGRI 
(see Figure|2j, and best fit parameters of the average spectrum of AX J1841. 0—0536 observed with IBIS/ISGRI (see Figure|3j. F is the powerlaw photon 
index, Ec is the cutoff energy, Bp is the e-folding energy. 



Source 


Fit model 


A^H (10^2 cm-2) 


r 


Be, Bp (keV) 


xl (d.o.f) 


JEM-X & IBIS/ISGRI 


IGR J 16479-45 14 
XTE J 1739-302 
IGR J 17544-26 19 


cutoff-powerlaw 
cutoff-powerlaw 
powerlaw with high-energy cutoff 


23 + 1'' 
8.6+" 

— 3.6 
2+* 


1 2+"-'* 
1 6+°-2 

^ ^+0.3 
^■^ -0.1 


Be = 25.ltig 
Be = 26.5+6-5 

Be = 19.7+i-^;Bp = 8.1+J'° 


0.90 (23) 
1.12 (45) 
1.03 (27) 


IBIS/ISGRI 


AXJ1841.0-0536 


powerlaw 




2 50 +"■1" 




1.20(10) 



the approximated formula obtained bv lEggletonI ( Il983h . is the fol- 
lowing: 



0.6g2/3 +ln(l + gi/3) ^ ^ 



tic0, with a circular orbit assumed. Figure |5] reports the results 
obtained for _Rl. Roche Lobe overflow (RLO) is not expected in 
this system because RLO would imply a much higher mass trans- 
fer to the compact object and thus a much higher X— ray lumi- 
nosity. Assuming that the compact object is a neutron star with 
A/x — 1.4 Mq, we make the hypothesis that the primary star has 



2 Mp is the mass of the primary star (i.e. the supergiant star), and A/x is 
where a is the orbital separation, q = Mp/M^ is the mass ra- the mass of the compact object. 
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Figure 5. Rl for different values of q = Mp/M^ (see equation[T). 
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Figure 6. Comparison between observed (solid line) and calculated (dashed 
line) distributions of the flare luminosities of IGR J16479— 4514. 



mass Mp = 31 Mq and radius Rp — 19 Rq. This set of parame- 
ters is in agreement with the expected values for a 08.5I star (e.i 
iMartins. Schaerer & Hillie3l2005l : l\^cca. Garmanv & ShullllT996 
and excludes the Roche Lobe overflow. However, we point out that 
other parameters for Afx, Mp, and Rp are allowed. 

We a pplied our clumpy wind model in its spherical con- 
figuration l lDuccietal.ll2009f) to the IBIS/ISGRI observations of 
IGR J16479— 4514, for which we have a significant number of 
flares (Na = 80). We will apply this model to the other SFXTs in 
a future work, when the number of their flares will be enough. Our 
model assumes that a fraction / of the mass lost by the supergiant 
wind is in form of clumps (/ = Mci/Mtot)- The clumps are as- 
sumed to follow a power-law mass distribution p(AIci) oc M^^^ in 
the mass range Afa— Afb, a power-law radiu s distribution A*' oc 
and a /3-velocity law (see lDucci et al J2009l for details). 

We compared the observed and calculated number of flares, 
their luminosity distribution, the X— ray luminosity outside flares 
and the average time duration of flares. 

For each of the 80 flares found, we derived the peak luminosity 
in the energy range 1 — 100 keV by means of the spectral param- 
eters found by fitting simultaneously the IBIS/ISGRI and JEM-X 
data (see TableQ. The lumi nosity (1 — 100 ke V) outside flares is 
Lout-flaros < 10^^ erg s'^ jSgucra ct alj|2008h . 

To avoid the selection effect given by the fact that the 5a de- 
tection threshold varies with the position of the source in the field of 
view, we have considered only those flares with luminosities greater 
than Liim ~ 4 x 10"^® erg s~^, where we have calculated Lum from 
the average 5a detection threshold count— rate at an offset angle 
from the pointing position of 6 = 15°. 

Figure[6] shows the comparison between the observed and cal- 
culated distributions of flare luminosities which results in the best 
agreement. The wind parameters obtained are reported in Table 
m with M = 10"^ Mq yr~\ v^c = 1800 km s-\ /3 = 1, 
Ma = 5 X 10^^ g, Mb = 5 X 10^1 g, C = 1.1, / = 0.5, 7 = -4, 
and, since Teff = 34000 K, the multiplier parameters adopted from 
Shimada et all (1994) are: k = 0.375, a = 0.522, S = 0.099 (see 



Ducci et al I2OO9 for details). With this set of wind parameters we 



4 DISCUSSION 

The wind parameters obtained for IGR 116479—4514 in section 
13.11 are in agreement with thos e of the source s previously stud- 
ied (see Ducci et al. 2009; Rom ano et al.ll2O10l) . with the excep- 
tion of / (/ = 0.75 for the other persistent HMXBs and SFXTs). 
Moreover, the mass loss rate found for IGR 116479—4514 (M = 
10~^ Mq yr~^) is lower if compared to the typical mass loss 
rate from a 08.51 star, which is of the order of ~ 4 x 
10"" Mq yr~^ (see Table [5). This difference could be due to 
the fact that the mass loss rates derived from homogeneous-wind 
model measurements with the Ha method are overes timated by a 
factor 2-10 if the wind is clumpy (see e.g. iLepine & M offat 200^; 
iHamann. Fe ldmeier & Oskinova l'2008). Another possibility is that 
the X-ray flaring behaviour in IGR J16479— 4514 is not totally 
due to the accretion of clumps, but other mechanisms could be at 
work, like, e.g. the centrifugal inhibition of accretion, the formation 
of transient accretion disks, or the Rayleigh— Taylor instability. The 
role of these mechanisms in SFXTs will be treated in this section. 



4.1 X-ray photoionization 

4.1.1 Direct accretion 

Here we discuss the effects of X-ray photoionization on the mass 
transfer onto a neutron star in a close binary system, modifying the 
analytic model bv lHo & AronsI ( 1 19871) and applying it to SFXTs. 

In HMXBs, the wind reaching the compact object is overi- 
onized by the X-ray pho tons produced by the compact object 
l lHatchett & McCravlll977h . This high ionization alters the dynam- 
ics of the line-driven stellar wind of the primary: the wind becomes 
highly ionized and does not interact anymore with the UV pho- 
tons emitted by the primary, which drive the wind. Hence, the wind 
velocity and density around the neutron star is different from what 
expected neglecting the X— ray photoionization. The accretion onto 
the compact object is modified accordingly since it depends onto 
the wind velocity. 

For a source with a X— ray luminosity Lx, the ionization state 
5 of the wind at a radius R^ is defined by: 



are able also to reproduce the observed X— ray luminosity outside 
bright flares and the average flare duration. 



(2) 
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Figure 7. Schematic representation of a Stromgren surface produced by 
a neutron star orbiting around a supergiant. The labels indi cates the dis- 
tance o f this surface which respect the two stars, accordingly to lHo & Aronj 
h987l) . 
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Figure 8. Schematic representation of the wind velocity, described by 
equation (5), and the hmit-velocity (equation 112) . We assumed v^o = 



1800 km s" 



1, -Rd 



19 Rp, and M„ = 31 Mp,, M, 



1.4 Mr: 



where n{r^) is the particle number density at a distance = 
D — from the supergiant (see Figure |7}. The parameter ^ 
defined in equation ^ determines the thermal and ionic state 
of the gas, assuming optically thin gas in thermal balance (see 
iTarter. Tucker & Salpeteill969h . 

Several attempts have been made to include the ef- 
fects of X— ray io n ization in the accret i on of HMXBs (e.g. 
Ho& AronsI 1 19871; iMacGregor & Vitelld 1 19821 : .Blondin et al., 



1990l : IStevens & KallmanHl99(]|:IStevenslll99l[). H ere we adopt the 



analytic model developed bv^ Ho & AronsI ( Il987[) . They assumed a 
spherically symmetric wind, ionized by the X— rays from the com- 
pact object. They assumed that the radiation line force is turned off 
at ^cr ~ 10* erg cm s~^. Thus, the wind follows the standard 13- 
velocity law (with 13 — 0.5) up to a distance Rcr from the compact 
object such that ix/[7i(i?cr)iicr] = Ccr, i.e. where the wind is suf- 
ficiently ionized to become transparent to the UV photons. Inside 
this sphere with the neu tron star in the centre, the radiation force 
jKudritzki & PulslllOOOl) is turned off and the wind velocity is as- 
sumed constant. The wind velocity will be lower in the vicinity of 
the compact object (with respect to the non-ionized case), leading 
to an enhancement of the mass accretion rate. With these assump- 
tions, IHo & AronsI ( Il987l) developed a model to describe the ac- 
cretion of the wind, taking into account the feedback effect of the 
X— rays which ionize the wind, and consequently controls the mass 
transfer onto the compact object. 

We propose the following important changes to this model: 

(i) a generic /3-velocity law, with /3 not fixed; 

(ii) we consider the orbital velocity of the neutron star in the cal- 
culations, which cannot be neglected i n close binary s ystems such 
as IGR J16479-4514 (Poa = 3.32 d, |jain et al.ll2009h : 

(iii) we modified the equations developed by IHo & AronsI 
l ll987l) to take into account the possibility that the mass loss rate 
towards the neutron star is reduced because of the high ionization 
state of the wind, which reduces the radiative acceleration given 
by the absorption and re-emission of UV photons (emitted by the 
supergiant) in the resonance lines of ions forming the wind. Assum- 
ing for the wind particles a velocity distribution centered on i;w (r^) 
(equation [3j, a part of them may decelerate enough not to be able 
to reach the neutron star; 



(iv ) a force cutoff value ^cr = 3 x 10^ erg cm s ^ ( IStevensI 
fT99lh . 

We assumed the standard win d velocity law obtained from th e 
radiation line-driven mechanism of lCastor. Abbott & Klein! J 19751) : 



(3) 



where is the distance from the primary (see Figure|7j, and Voo is 
the terminal velocity. 

In wind-fed systems, the mass accretion rate is defined as the 
flux of matter passing through a circular area with radius i?a, and 
is written as: 



M 



4:■KD^V^ 



-Vrcl 



(4) 



where M is the wind mass loss rate from the supergiant, the factor 
contained in square brackets is the stella r wind flux at a distance D 
from the primary, v^ci ~ -^/t^w + f orb relative velocity be- 

tween the wind and the neutron star, and iia is the accretion radius, 
defined as: 
_ 2GM. 

where Mx — 1.4 Mq is the mass of the neutron star and Cs is the 
sound velocity. Hence, the accretion luminosity, given by equations 
lUlElllll is: 



GMy, 



Ma, 



1 



(6) 



Px i?x DH^{r^){vl, + af/^ 

The radius riim where the sum of the radiation force g^, due 
to scattering of continuum photons by free electrons, and the gravi- 
tational forces of the primary and neutron star acting on a parcel of 
the wind is zero (see Figure[8]l can be derived as follows: 

GAfx 



where Fc, defined as 



(7) 



(8) 



AncGMp 

is obtained from the formula of the force due to electron scattering: 
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GM, 



J- c 



(9) 



where CTo is the opacity for electron scattering, assumed to 
be equal to ~ 0.3 cm^ g^^ and constant, as suggested by 
iLamers &CassinelliHl999l) . 

When the radiation force is turned off, a parcel of gas is only 
subject to the gravitational force of the two stars and to the force 
due to electron scattering. Assuming that this parcel of gas reaches 
the distance rum with null velocity, the initial velocity can be ob- 
tained from the momentum equation: 



dv 
dr 



GAL 



GM^ 
(D-r) 



(10) 



For a particle starting at a distance r, with velocity v, the solution 
obtained from equation l IlOl l is: 



GAU 



10 

Thus 



-2GMp(l-re 

1 



2GMx 



D-r 



(D 



1/2 



rdr'.(ll) 



(12) 



The wind particles have a temperature of ~ IC K 
dLamers & Cassinellill999ll . Assuming for the particles velocity, at 
a distance r from the primary, a gaussian distribution centered on 
Voc{i — Rp/r)^ with CTv ~ 10^ cm s^^, it is possible to calculate 
the density probability to have particles with v 



air — rc 



27r(Tv 



dv 



(13) 



where u„ is given by equation JS}. 

We now propose a new version of the self- consistent steady 
state equation developed by ' Ho & Aron j imt . improved with 
the considerations described above. The accretion luminosity La 
from equations l[3]l, 1(6), and l ll3t is then: 



(GMx)' 



(14) 



which we call, following the nomenclature o f lHo& Aronsl l ll987h . 
the accretion equation. Then we obtain the feedback equation from 
equations (O, ig}, lUll, (THi: 



Lh{r{) = 



MS{r^){D-r^) 



4nr'j^mpVc. 



\l-Rp/r,y 



(15) 



where /i is the mean atomic weight. The steady state solution is 
given by: Lii{r^) — L\^{r^), thus, by equating equations il4\ and 
il5i , it is possible to obtain the X— ray luminosity Lx and the cor- 
responding Stromgren radius ra- 

We assume for the force cutoff value ^cr ~ 3 x 
10^ erg cm s~^, which is two orders of magnitude smaller than 
the value considere d bvlHo & Aronsi ( Il987 ). but in agreement with 
the calculations of StevensI jl99lh . who found that for this value the 
wind material is already basically completely ionized. 

We calculate the X— ray luminosity and rcr of 
IGR J 16479— 4514 from the steady state solution obtained 
from equations J14t and il5\ . assuming the parameters reported in 
Table|5] We obtain an expected luminosity of Lx ~ 10^^ erg s~^, 
reached by the neutron star when rcr ~ 1-33 x 10^^ cm. This 
luminosity level is in agreement with the peak flare luminosities 



Table S. System Parameters for IGR J16479-4514 and IGR J17544-2619. 



Parameters 



Sources 

IGRJ16479-4514 IGR J17544-2619 



Spectral type 


08.5I 


091 


forb (d) 


3.3194 ± 0.0010*^ 


4.926 ±0.0018 


Mp (Mq) 


31 " 


25 - 28 s 


Rp (Rq) 


19" 


< 23 if e = s 






~ 12.7 if e 0.4 s 


Lp (L©) 




Si 4.6 X 10^ 


M(Mq yr-1) 


~ 4 X 10-'^ " 


2.4 X 10"*^ " 


Mx (Mq) 


1.4 t 


1.4f 


Voa (km s"'^) 


1800 


1800 


/3 


1 


1 



^ iRahouiet al]j2008l) 
IjainetalJiioO^ 

obtained assuming no Roche-Lobe overflow (see section [3Tl 

IVacca et alj h996l) : lMartins et al.l f200l 
" Calculated with the Vink formula jVink, de Koter & Lamersl2000l) 
^ Assumed 
s IClarketalj|2009h 

Table 6. Results for IGR J17544-2619. 





Periastron 


Apastron 


M(Mq yr-1) 


2.4 X 10^** 


2.4 X 10-6 


rcr (cm) 


~ 0.9 X 10^2 


~ 3.3 X 10^2 


Lx (erg s-l) 


> 1037 


~ 2 X 10^5 


tiroi (km s-i) 


~ 230 


~ 1400 


Rm (cm) 


lOio 


4.5 X 



observed. However, IGR J 16479— 4514 h as an out-of-flare lumi- 
nosity of the order of lO^'' - 10^^ erg s"^ jSidoli et alj2008h . This 
luminosity level can be obtained possibly invoking the presence of 
the centrifugal inhibition of accretion. 

For IGR J17544— 26 19 we assumed the system parameters 
suggested by IClark et al.l ( 12009.) , with eccentricity e = 0.4. We 
found that the X-ray photoionization acts to increase the differ- 
ence in luminosity when the neutron star is at periastron (Lx 



10 



•^^ erg s ^) and apastron (Lx 



10^^ erg s-\ see Table[6l(. In 
fact, at periastron the wind is highly ionized by the X— ray source 
(rcr — 0.9 X 10^^ cm), leading to a low wind velocity in the vicin- 
ity of the neutron star (iirci — 230 km s"^) and accordingly to a 
higher X-ray luminosity (since Lx oc ). At apastron the ioniza- 
tion is lower (rcr — 3.3 x 10^^ cm), therefore its effect on the wind 
velocity is small (Urd — 1400 km s"^). 



4.1.2 Centrifugal inhibition of accretion 

In this section we apply for the first time t he centrifugal inhib ition 
(c.i.) of accretion to the modified model of IHo & Aronsi ( Il987h de- 
scribed in section [4.1.11 

The magnetospheric radius Rm is defined as the radius where 
the magnetic field pressure i?^(7?m)/87r equals the ram pressure 

We use the definition ob- 



of the ac creting plasma p(Rm)v^ (Rm 
tained bv lDavidson & Ostrikeilll973l) : 



i?m(re) 



BqRx 



4Afaccr(GMx)l/2 



2/7 



(16) 
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where Bo is the surface magnetic field. When 7?m > _Rco (where 
_Rco is the corotation radius), the accretion flow is halted at the 
magnetospheric boundary, which behaves like a closed barrier 
dlUari onov & Sunvaev. .1975i) . The expected X— ray luminosity in 
this regime is: 

Lm = ^^Maccr . (17) 

The position of the magnetospheric radius depends also on the ef- 
fect of the X— ray ionization on the accreting matter. Here we ob- 
tain the formula which gives the steady position of the magneto- 
spheric radius and the related position of the Stromgren radius for 
which the accretion luminosity is equal to the feedback luminosity. 
The accretion equation il4l can be written as follows: 



(GMx) 



MS{r^) 



Rn 



2-13/2 



(18) 



We obtain Rm as a function of rcr by equating Lb of the feedback 
equation il5\ to the La.c.i. of the accretion equation JlSt : 



-Rm,c.i.(r^) 



(GMx)^47r^mp 



(19) 



Hence, if the centrifugal inhibition of accretion is at work in an 
X— ray binary system where the X— ray ionization cannot be ne- 
glected, the steady values for Rm and rcr can be derived by equat- 
ing equations il6l and l |19l l: 



(20) 



In Table |7] we report the steady state solutions obtained for 
IGR J16479— 4514 (assuming the system parameters reported in 
Table Is} in the case of direct accretion and in the case of centrifu- 
gal inhibition of accretion. Assuming a transition from direct accre- 
tion to the centrifugal inhibition of accretion for a mass-loss rate of 
M = 4 X 10~® Mq yr~^, we obtain steady state solutions in agree- 
ment with observations for a magnetic field of Bo ~ 10^^ Gauss 
and a spin period of ~ Is. When the neutron star enters the state 
of centrifugal inhibition of accretion, the X— ray luminosity is re- 
duced, thus also the photoionization is lower; this leads to a shift of 
the magnetospheric radius to a higher value. The opposite transition 
(from centrifugal inhibition to direct accretion) needs A/ greater 
than 4 x 10~^ Mq yr~^ (e.g. the accretion of a dense clump). This 
behaviour allows longer durations for the low-luminosity state. 
When i?ni ~ i?co, t he oblate spheroid shape of the m agneto- 
spheric boundary (see Ijetzer. Strassle & Straumann|[l998h results 
in the intermediate luminosity state observed (~ 10'^^ erg s~^) be- 
cause simultaneously Rm < Rco in the mag netic polar region and 
Rm ^ Rco in the magnetic e quatorial region dCampana et alj20()ll : 
IPema. Bozzo & Stellall2006h . 

4.2 Formation of an accretion disk 

In the previous section we discussed the effect of X— ray photoion- 
ization in a close X— ray binary system like IGR J16479— 4514, 
where a compact object accretes matter from a wind with a low 
relative velocity v^ei ~ 460 km s~^. In section [4.1. 11 we have 
excluded for this system the possibility of the formation of an 
accretion disk due to Roche Lobe overflow. However, it is also 
possible that an accretion disk form s from the capture of angu- 



lar momentum from a slow wind (see Illarionov & Sunvaev|[T973 : 



Table 7. Model results for IGR J16479-4514. We calculate the magne- 
tospheric radius Rm with equation )16t . assuming a magnetic field of 
Bq = 10^^ Gauss. For IGR J16479-4514, the corotation radius is equal 
to the magnetospheric radius when M = 4 X 10~® yr~^ for a spin 
period of the neutron star of ~ 1 s. 





Direct accretion 

{Rm < Rco)'- 


Centrifugal inhibition 
of accretion (-Rm > Rco)'- 


M(Mq yr-1) 


> 4 X 10-" 


4 X lO-'' 


rcr (cm) 


~ 1.33 X 10^2 


~ 2 X 10^2 


Lx (erg s~i) 


> 1037 


> 8 X 10^4 


i)rci (kms-i) 


~ 460 


~ 760 


Rm (cm) 


1.6 X 10* 


1.7 X 10* 




1x10" 2x10" 3x10" 4x10" 5x10" 
time (s) [start time = 53612.164 MJD] 



Figure 9. Lightcurve of the observation of XTE J1739— 302, in the energy 
range 18 — 60 keV (IBIS/ISGRI). Arrows indicates the peaks of luminosity. 



formula [31] in IWand[T98lh . therefore the presence of a disk in 
IGR J 16479— 4514 caimot be ruled out. The formation of a tran- 
sient disk from the mass and angular momentum capture from an 
asymmetr ic ste l lar wi nd can lead to a flaring activity, as proposed 
by Taam e t alj il988|) to explain the recurrent flares observed in 
EXO 2030-1-375 jParmaretal.1 1 1989h. Recently, the formatiori of 
transient disks has been proposed by Krevkenbohm et al] l l2008l) to 
explain part of the fla ring behaviour of t he persistent HMXB Vela 
X— I. In the model of iTaam et alj ( Il988[) . the density and velocity 
inhomogeneities in the wind lead to an instability in the accretion 
flow. Because of this instability, the interaction of the accretion flow 
with the shock fronts of the accretion wake leads to the formation 
and dissipation of transient accretion disks (see e.g. Edgar 2004 and 
references therein for a recent review). In their model. ITaam et all 
(1988) gave an approximative formula for the time scale for the 
duration of flares: 



6GMx 



(21) 



IShapiro & Lightmanlll976l : IWand f I98ll) . The existence of a disk the duration of flares of 
in IGR JI6479— 4514 requires that the relative velocity between 
the neutron star and the wind is lower than ~ 500 km s"^ (see 



Assuming relative velocities reported in Table|7](in the case of di- 
rect accretion), we found for IGR JI6479— 4514 a time scale for 



10'' 



10 s, in agreement with the 



observed flare durations. The flares produced by the formation of 
transient accretion disks can be distinguished from those produced 
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Figure 10. Lightcurve of an observation of IGR J16479— 4514, in the en- 
ergy range 18 - 60 keV (IBIS/ISGRI). 



by the simple clump accretion by measuring the time derivative of 
the pulse period, which, in the first case, is expected to change sign 
for each flare fsee lTaam et alJ[l988h . 

The SFXT XTE J1739— 302 has shown a quasi-periodic flar- 
ing behaviour on 53612.164 MJD (see Figure|9]l, which could be 
ascribed to the mechanism of formation and dissipation of tran- 
sient accretion disks described above. Adopting equation | |2U . we 
obtain, from the measure of the flares durations (r ^ 5000 s), a 
relative velocity of t^rci ~ 600 km s"^. This value is in agree- 
ment with the Urci expected from an X— ray binary system with a 
small orbital period (« 3 — 7 d), or with a higher orbital period 
and large eccentricity, where the effect of X— ray photoionization 
reduces significantly the wind velocity. This supports the hypothe- 
sis of Smith et al. ( 2006 ) that the fast outbursts of XTE J1739-302 
could be due to an instability of an accretion disk. 



4.3 Accretion via Rayleigli— Taylor instability 

The shape of the flare of IGR J 16479— 4514 we observed with 
INTEGRAL and reported in Figure [TO] i.e. fast rise and exponen- 
tial decay, could be explaine d with the magne tospheric instabil- 
ity mechanism proposed by iLamb et al.l jl977h to explain type II 
X— ray bursts. All equations in this section are taken from these au- 
thors. In the framework of a spherically symmetric accretion flow 
onto a neutron star, and under particular conditions of X— ray lu- 
minosity and of the temperature of the accreting matter, the mag- 
netospheric surface behaves as a gate which controls the flow to- 
wards the stellar surface. When the gate is closed, a magneto- 
spheric cavity is formed, and a reservoir of matter is accumu- 
lated on the top of the neutron star magnetosphere, leading to low 
X— ray luminosities. When this matter has cooled enough (because 
of electron— ion bremsstrahlung, which is the dominant cooling 
mechanism when the gate is closed), the plasma enters the magne- 
tosphere via Ray leigh— Taylor instability, and the accretion onto the 
neutron star leads to a flaring behaviour. iLamb et alj tl977i) found 
a critical luminosity 



2x10'' 



M 

Me 



1/2 



i?6-'/'erg s-i(22) 



where Tr is the temperature of radiation, Tc ~ 10^ K, /130 is the 
stellar magnetic moment in units of 10^" Gauss cm^, M is the mass 
of the neutron star and Ra is the radius of the neutron star in units of 
lO'' cm. If the flare luminosity is greater than the critical luminosity 
(equationl22t. the Compton cooling dominates the bremsstrahlung 
cooling at _Rni, thus the m agnetosphere gate is open for longer time 
and the flare is prolonged jLamb et"al]|l977h . 

The flare reported in Figure[lO]has a peak luminosity of ~ 3 x 
10'^^ erg s~^, which is greater than the critical luminosity Lcrit ~ 
2.4 X lO^'' erg s~^ (if < 10^ K). 

Assuming Tr ~ 10* K, we found that, because of Compton 
processes, beyond the radius ~ 2.2 x 10^" cm at which the free- 
fall temperature equals the temperature of radiation, the X-rays 
emitted by the ne utron star heat the plasm^ (see formula [31] in 
lLambetalJll977l) . 

If the heating time scale is less than the flow time scalfl the 
flow is choked at i?m, and we expect for the flare a time duration 
given by: 



125 



-3/2 



M 

Mq 



(24) 



which results into Afmax ~ 175 s assuming Tr = 10* K and M = 
1.4AfQ. This duration is not in agreement with that observed in 
Figure[TO] However, if we assume Tr ~ 2 x 10^ K, the heating time 
scale is less than the flow timescale at distances 3> r^. The condi- 
tion l l23t is no longer valid (Lx < Lchokc ~ 4 x 10^*^ erg s~^), 
thus the flow is choked at a radius r'^ > fc. In this case the expected 
burst duration is given by: 



At' 



210L7, 



10* if 



M 



s ^ 4900 s 



(25) 



(assuming = 2 x 10^ K, M = I.4M0 and Lx = 3 x 
lO''^ erg s^^). The assumption Tr = 2 x lO'^ K is in agreement 
wit h the temperature of radiation observed in some SFXTs (see 
e.g.lSidohetalJl2007l) . In the framework described bv lLamb et al.l 
l ll977 ). the exponential decay in this kind of flares is described by 
the law L{t) oc t^"^"/^, with a > 3/2. From the fit of the ex- 
ponential decay of the flare observed by INTEGRAL (Figure [TOt 
we found a ~ 1.5 5, thus in agreement with the value predicted by 
lLambetalJ ( ll977h . 

We point out that if Lx > Lcrit, Tr < Tc at Rm, and if 
the accretion radius is lower t han r^, the X-ray sou rce is persis- 
tent with high luminosity (see lElsner & Laniblll977l) . In our case 
7?a ~ 1.8 X 10" cm and rc ~ 2.2 x 10^° cm (if = 10* K) 
or rc ~ 1.1 X 10" cm (if = 2 x lO'^ K). Thus, what makes 
IGR J 16479-45 14 (and probably other SFXTs) intermittent with 



The material accreting onto a neutron star is decelerated by the magnetic 
field, hence it is shock-heated to a temperature Ts ~ 3/16Tff (Rm), where 
Tff (r) oc 1/r is the proton free-fall temperature. Since the temperature of 
radiation is roughly constant for every r, below a radius rc the gas is cooled 
by Compton interaction (Tr < Tgas); outside rc we have T > Tgas, thus 
the gas is heated by Compton processes. 

* The condition for which the heating time scale is less than the flow time 
scale is expressed by: 



Lx > Lc 



2 X 10^'^ ( ( — ] erg s'^ . (23) 



This condition is respected if T = 10* K, M = 1.4Mq, = 3 x 
10^^ erg s-^ 
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respect to the other persistent HMXBs could be the X— ray pho- 
toionization which reduces the wind velocity, increasing the accre- 
tion radius with respect to rc. 



5 CONCLUSIONS 

Up to now three accretion mechanisms have been proposed to ex- 
plain the X— ray flares of SFXTs. 

Here we propose other mechanisms to explain the observed 
X-ray behaviour of IGR J16479-4514, IGR J17544-2619 and 
XTE J1739— 302: the effect of X— ray photoionization onto accre- 
tion (both in direct accret ion and centrifugal i nhibition of accretion) 
in the framework of the IHo & AronsI ( Il987b model; the accretion 
disk instab ility of Taam et ^T i 1988h . and the Rayleigh-Taylor in- 
stability of lLambetal.l ( ll977h . 

We have shown that the X— ray photoionization reduces the 
mass loss rate and the wind velocity along the trajectory compact 
object— OB supergiant with respect to the undisturbed case (section 
|4.1| l. Their simultaneous reduction leads to X-ray luminosities in 
agreement with those observed during the flaring activity; more- 
over, a lower allows the formation of transient accretion disks 
from the capture of angular momentum, able to reproduce some 
kind of quasi— periodic recurrent flares observed in SFXTs. 

We found in our analysis of INTEGRAL data that some flares 
show a peculiar shape, characterized by a sharp rise and an ex- 
ponential decay (Figure [TOt. W e found that t his sh ape is well ex- 
plained assuming the model of lUamb et iD jl977h . based on the 
Rayleigh-Taylor instability, although this model was proposed to 
try to explain type II bursts. 

Therefore, we conclude that in SFXTs with large orbital peri- 
ods (Porb > 15 d; e.g.: IGR J18483-0311, SAX J1818.6- 1703, 
IGR Jl 1215-5952) the effects of X— ray photoionization onto 
accretion mechanism can be neglected. In SFXTs with smaller 
orbital periods (Porb < 15 d; e.g.: IGR J16479-4514, 
IGR J17544-2619) we expect that part of the X-ray variability 
observed is due to the X-ray photoionization and consequently 
to the accretion disk instability. We are still not able to establish 
which mechanism dominates in SFXTs with smaller orbital peri- 
ods. Moreover, we point out that different accretion mechanisms 
could be at work in a single SFXT, in the case of orbits with high 
eccentricities. 
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